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Abstract 

To better understand and manage complex systems like ecosystems it is critical to know the relative con- 
tribution of the system components to the system function. Ecologists and social scientists have described 
a diversity of ways that individuals can be important; This paper makes two key contributions to this re- 
search area. First, it shows that throughflow (Tj), the total energy or matter entering or exiting a system 
component, is a global indicator of the relative contribution of the component to the whole system activ- 
ity. Its global because it includes the direct and indirect exchanges among community members. Further, 
throughflow is a special case of Hubbell status or centrality as defined in social science. This recognition ef- 
fectively joins the concepts, enabling ecologists to use and build on the broader centrality research in network 
science. Second, I characterize the distribution of throughflow in 45 empirically-based trophic ecosystem 
models. Consistent with theoretical expectations, this analysis shows that a small fraction of the system 
components are responsible for the majority of the system activity. In 73% of the ecosystem models, 20% 
or less of the nodes generate 80% or more of the total system throughflow. Four or fewer nodes are required 
to account for 50% of the total system activity and are thus defined as community dominants. 122 of the 
130 dominant nodes in the 45 ecosystem models could be classified as primary producers, dead organic 
matter, or bacteria. Thus, throughflow centrality indicates the rank power of the ecosystems components 
and shows the concentration of power in the primary production and decomposition cycle. Although these 
results are specific to ecosystems, these techniques build on flow analysis based on economic input-output 
analysis. Therefore these results should be useful for ecosystem ecology, industrial ecology, the study of 
urban metabolism, as well as other domains using input-output analysis. 

Keywords: input-output analysis, food web, trophic dynamics, social network analysis, ecological network 
analysis, materials flow analysis, foundational species 



1. Introduction 

Identifying functionally important actors is a critical step in understanding and managing complex 
systems, whether it is a fortune 500 company or an ecosystem. For example, Ibarra (1993) showed that an 
employee's power to affect administrative innovation within an advertising agency was in part determined by 
their positional importance within the organization. In ecological systems, knowing the relative functional 
importance of species or groups of species is essential for conservation biology, ecosystem management, 
and understanding the consequences of biodiversity loss (Walker, 1992; Lawton, 1994; Hooper et al., 2005; 
Jordan et al., 2006; Saavedra et al., 2011). 

Ecologists have several ways of classifying the relative importance of community members. Whittaker 
(1965) introduced rank-abundance curves to describe the community richness and indicate the relative 
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importance of the species, assuming that community importance was proportional to abundance. He also 
presented an alternative rank-productivity curve that indicated the species importance based on their net 
productivity. Subsequent ecological concepts have built on this. Keystone species (Paine, 1966; Power et al., 
1996) are species whose importance to the community are disproportionate to their biomass, like the sea 
otter in Pacific kelp forests. Ecological engineers (Lawton, 1994; Jones et al., 1994) are species whose actions 
create whole new habitats, such as beavers that transform terrestrial environments into slow moving aquatic 
environments. Dayton (1972) introduced the more general term foundational species for fundamentally 
important species of many types (Ellison et al., 2005). Part of the challenge and the reason for multiple 
concepts, is that there are a diversity of ways in which a species may be important and contribute to a 
community or ecosystem. 

Faced with the analogous problem of identifying important members of human communities, social 
scientists developed the centrality concept (see Wasscrman and Faust, 1994). Centrality embodies the 
intuition that some community members are more important, have more power, or are more central to 
community function. Centrality was developed in the context of network models of communities in which 
individuals are represented as nodes of a graph and the graph edges signify a specific relationship between 
two individuals such as friendship or co-authorship (Fig. la). The relationship may or may not be directed. 
Degree centrality is the number of immediately adjacent neighbors on the graph, and it assumes that more 
connected nodes are more central. It is quantified as the number of edges incident to the node. In the 
example graph, node 3 has a degree of 7 (note the separate directed pathways from 3 to 4 and from 4 to 
3 shown as a two headed arrow). Fig. lb shows the distribution of node degrees in the community which 
indicates that node 3 is the most central from this local neighborhood perspective. 

Some scientists have suggested that the local neighborhood is insufficient to determine the node's cen- 
trality for some applications, especially exchange networks (Hubbcll, 1965; Bonacich, 1972; Estrada, 2010). 
Instead, a node's importance may be increased because one or more of its neighbors are important. Network 
models can capture this increased neighborhood size by defining a walk as a sequence of edges traveled 
from one node to another, and walk length (m) is the number of edges crossed. In the example network, 
there is a walk from 6 to 2 of length to = 3 by following 6 — >■ 4 3 — >■ 2. This enables us to consider 
the neighborhood to steps aways (Estrada, 2010). Fig. Ic shows the eigenvector centrality (Bonacich, 1972, 
1987) for the example network which identifies the equilibrium number of paths passing through each node 
as to ^ oo. In this sense it is a global centrality measure because it is a "summary of a node's participation 
in the walk structure of the network" (Borgatti, 2005) and captures the importance of indirect as well as 
direct interactions (Borgatti, 2005; Scotti et al., 2007). 

Degree and eigenvector are only two examples of centrality indicators. Many centrality measures have 
been developed and applied in the literature for complex systems modeled as networks (Wasscrman and Faust, 
1994; Koschiitzki ct al., 2005). The centrality measures tend to be correlated (Newman, 2006; Jordan et al., 
2007; Valente ct al., 2008), but the differences can be informative (Estrada and Bodin, 2008; Baranyi ct al., 
2011). Borgatti and Everett (2006) provide a classification of centrality indices and shows how and why 
different measures are useful for different applications. 

Ecologists have applied the centrality concept in several ways. For example, landscape ecologists have 
used centrality to assess the connectivity of habitat patches, how this connectivity effects organism move- 
ment, and how habitat loss changes the connectivity (Estrada and Bodin, 2008; Bodin and Saura, 2010; 
Baranyi et al., 2011). Community and ecosystem ecologists have developed and used centrality measures to 
study how organisms influence each other in transaction networks (Jordan ct al., 2003; Allesina and Pascual, 
2009; Fann and Borrett, 2012). Jordan et al. (2006) argue that mesoscale measures, between local and global 
centralities, are most useful for ecosystem studies because the impact of indirect effects tend to decay rapidly 
as they radiate through the system. Recent work used centrality indicators to determine important species 
in communities of mutualists (Martin Gonzalez ct al., 2010; Sazima et al., 2010). Collectively, this work 
shows how a range of centrality indicators can be useful for addressing ecological questions. 

Here, I identify a new centrality indicator for ecology, termed throughfiow centrality Tj. I first recognize 
that the throughfiow measure ecosystems ecologists have long calculated (Patten et al., 1976; Finn, 1976; 
Ulanowicz, 1986) is a global measure of node importance in generating the total system activity. Further, 
I show that this is a special case of Hubbell's status index centrality (Hubbell, 1965). I then apply this 



63 measure to 45 trophic ecosystem models drawn from the hterature to test two hypotheses regarding ecosys- 

64 tem organization. The first hypothesis suggested by both Whittaker (1965) and Mills et al. (1993) is that 

65 communities are composed of a relatively few dominant species and larger group that are less central. The 

66 second hypothesis is that in ecosystems the dominant species/groups are expected to be comprised of pri- 

67 mary producers, decomposers like bacteria, and non-living groups included in ecosystem models like dead 

68 organic matter. This hypothesis stems from trophodynamic theory and energetic constraints of food chains 

69 (Lindeman, 1942; Odum, 1959; J0rgensen et al., 1999; Wilkinson, 2006) 

70 2. Theory — Throughflow is a Centrality Indicator 

71 A core claim of this paper is that the amount of energy-matter flowing through each node j in an 

72 ecosystem network — termed node throughflow (T, ) - is a global centrality indicator of the node's functional 

73 importance. In fact, this centrality measure is a special case of Hubbcll's (1965) status score. Further, this 

74 centrality indicator is more useful for ecologists and environmental scientists than the classic eigenvector 

75 centrality or the recently introduced environ centrality (Fann and Borrett, 2012) because (1) it is more 

76 intuitive to calculate, (2) it integrates the transient and equilibrium effects as flow crosses increasingly 

77 longer pathways, and (3) it captures the effects of environmental inputs (outputs) on the system flows. This 

78 section provides evidence to support these claims. 

79 2.1. Flow Analysis 

80 Flow analysis is a major branch of ecological network analysis (ENA) (Patten et al., 1976; Finn, 1976; 

81 Ulanowicz, 1986; Schramski ct al., 2011). It is an environmental application and development of Leontief's 

82 (1966) macroeconomic input-output analysis flrst imported to ecology by Hannon (1973). It traces the move- 

83 ment of energy-matter through the network of transactions in an ecosystem to characterize the organization 

84 and development of the system. 

85 2.1.1. Model Definition 

86 Flow analysis is applied to a network model of energy-matter exchanges. The system is modeled as a 

87 set of n compartments or nodes that represent species, species-complexes (i.e., trophic guilds or functional 

88 groups), or non-living components of the system in which energy-matter is stored. Nodes are connected by 

89 L observed fluxes, termed directed edges or links. This analysis requires an estimate of the energy-matter 

90 flowing from node j to i over a given period, F„xri = [/yli hi ~ 1,2,. (note the column to row 

91 orientation). This flux can be generated by any process such as feeding (like a food web), excretion, and 

92 death. As ecosystems are thermodynamically open, there must also be energy-matter inputs into the system 

93 z„xi = [zi], and output losses from the system yixn = [Vi]- In some applications, outputs are partitioned 

94 into respirations and exports to account for differences in energetic quality, but this is not necessary in this 

95 case. For other analyses, it is useful when the amount of energy-matter stored in each node (e.g., biomass) 

96 is also reported, x„xi — [xi] (Fath and Patten, 1999). The necessary model data can be summarized as 

97 X = {F,z,y,x}. 

98 To validly apply flow analysis, the network model must meet two analytical assumptions. First, the model 

99 must trace a single, thermodynamically conserved currency such as energy, carbon, or nitrogen. Second, 

100 the model must be at steady-state for many of the analyses. This means that the sum of the energy-matter 

101 flowing into a node equals that exiting the node such that its storage or biomass is not changing. Fath et al. 

102 (2007) offer further suggestions for better ecosystem network model construction. 

103 2.1.2. Throughflow 

104 Given this model, we can apply flow analysis. The technique has a dual approach. The input oriented 

105 analysis pulls the energy-matter from the boundary outputs and mathematically traces the pathways (a 

106 sequence of m edges) used to generate them all the way to the boundary inputs. In contrast, the output 

107 oriented analysis pushes inputs into the system and follows their paths through the system to their boundary 
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108 loss. This paper focuses on the output oriented analysis to support the centrality claims for brevity and 

109 clarity; the input perspective provides similar support. 

no The first analytical step is to calculate the node throughflows (T„xi — [Tj], j = 1, 2, . . . , n). Finn (1976) 

111 showed that the input and output oriented throughflows can be calculated from the initial model information 

112 as follows: 

n 

T^ = ^/,fc + z, (z = l,2,...,n), and (1) 

k=\ 
n 

^r^E^-^+y^ (j = l,2,...,n). (2) 

fc=i 

113 At steady state, [T/"] = [T^""'] — T and the amount of energy-matter stored in the node (xj) does not 

114 change through time. 

115 Finn (1976) argued that the sum of the node throughflows, called total system throughflow {TST — 

116 X]j=i ^i)' ^ measure of the activity or size of the ecosystem functioning. Ulanowicz and Puccia (1990) 

117 interpret Tj as the gross production of the compartment. Thus, Tj is the contribution of the j*'* node to 

118 the whole system functioning or productivity. It is in this sense that throughflow is a centrality measure 

119 indicating the relative importance or contribution of each node. 

120 Fig. 2 shows an example of rank ordered Tj for the Gulf of Maine ecosystem network (Link ct al., 2008). 

121 This shows the larger functional importance of phytoplankton, large and small copepods, detritus, bacteria 

122 in this system. This matches with the theoretical expectation that primary production and decomposition 

123 tend to be the critical components of ecosystem functioning (Wilkinson, 2006), but it also points to the 

124 importance of smaller consumers in the Gulf of Maine. Notice the similarity of this presentation to the 

125 rank-abundance and rank-productivity curves that Whittaker (1965) introduced to compare the relative 

126 importance of plants in a community. Like those original curves, Tj suggests that in this system there 

127 are a few dominant or more important species and a long tail of functionally less critical species (e.g., 

128 Pinnipeds, Beleen whales, and pelagic sharks). The application section considers the generality of both of 

129 these patterns. 

130 To facilitate comparisons between centrality measures, it is useful to consider the node throughflow 

131 scaled by the total system throughflow (Tj/TST) such that Tj/TST = 1. While the rank-ordering is 

132 preserved, rescaling in this way eliminates the units and differences in total magnitude between systems or 

133 other centrality measures. This focuses on intensive system differences while ignoring extensive differences 

134 present without the rescaling. Rescaling centrality measures is common, though it can introduce its own 

135 challenges (Ruhnau, 2000). 

136 2.1.3. Path Decomposition 

137 Path decomposition of throughflow lies at the core of ENA (Finn, 1976; Fath and Patten, 1999; Borrett ct al., 

138 2010), and shows why Tj is a global measure of functional importance. It partitions the flow of energy- 

139 matter from the input (output) over paths of increasing length (number of directed edges, m = 0, 1, 2, . . . , oo) 

140 within the system required to generate Tj . Recall that local centrality measures focus on the connections to a 

141 node's nearest neighbors or a restricted neighborhood, while more global measures consider the relationships 

142 between all nodes within the system. 

143 Path decomposition of flow starts by calculating the output oriented direct flow intensities G„xn = [9ij\ 

144 from node j to i. These intensities are defined as 

5., EE U,/T°-\ (3) 

145 Here, gij is the fraction of output throughflow at donor node j contributed to node i. The gij values are 

146 dimensionless and the column sums of G must lie between and 1 with at least one column less than 1 

147 because of thermodynamic constraints of the original model (J0rgensen et al., 1999). 
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148 The second step determines the output oriented integral flow intensities N = [n^j] as 

oo 

N = ^ G" (4) 

m=0 




Boundary Direct Indirect 



149 where Inxn = G° is the matrix muhiphcative identity and the elements of G™ are the fractions of boundary 

150 flow that travel from node j to i over all pathways of length m. As the power series must converge given 

151 our initial model definition, the exact values of N can be found using the identity N = (I — G)~^ . The riy 

152 elements represent the intensity of boundary input that passes from j to i over all pathways of all lengths. 

153 These values integrate the boundary, direct, and indirect flows. 

154 We can use N to recover T as follows: 



T = Nz. (6) 

155 This suggests that the path decomposition of throughflow shown in equation (5) is a true partition of the 

156 pathway history of energy-matter in the system at steady-state. 

157 The path decomposition in equation (5) shows how the throughflows are a global measure of centrality 

158 because the observed throughflows are generated by energy-matter moving over all pathways of all lengths 

159 such that the whole connected system is considered, not just a local neighborhood. Notice that the im- 

160 portance of longer pathways is naturally discounted as energy-matter is lost as it passes through nodes 

161 in the path. This discount or decay rate varies among ecosystems and model types (Borrett et al., 2010). 

162 Multiplication of the integral flow matrix by the boundary inputs to recover Tj (equation 6) illustrates how 

163 the node throughflows capture the potential effect of heterogeneous boundary inputs known to be a factor 

164 in ecosystems (Borrett and Freeze, 2011). 



165 2.2. Hubbell's Status Score 

166 Before Hannon (1973) applied Lcontief's (1965) economic input-output ideas to ecological systems, 

167 Hubbell (1965) applied the formalism to social systems. In doing so, he created a centrality measure that 

168 is known as Hubbell status or Hubbell centrality. Although Hubbell's initial model was different than the 

169 ecological one presented in section 2.1.1, the analytical mathematics is parallel to that shown for throughflow 

170 analysis. 

171 Hubbell (1965) started by modeling the interactions between individuals in a community using a weighted 

172 sociometric choice matrix W = [my],(i,j = l,...,n), where Wij can be positive or negative and indicates 

173 individual j 's indication of the strength of relationship between him or herself and individual i. The integral 

174 relationship strength among the community members propagated across the whole set of pathways R are 

175 then determined as 

R = I + + + W'"' + . . . , (7) 

176 where I„xn is again the matrix multiplicative identity and W" is the strength of relationship between any 

177 two community members over paths of length m. When the series converges, we can flnd R exactly as 

178 R=(I-W)-1. 

179 Building off of this analysis, Hubbell (1965) defined the status score S — [Si] of member i as 

S = R X E (8) 

180 where E„xi = [e^] are the system exogenous inputs. 

181 While the initial model was different, the throughflow equation (6) is identical in form to Hubbell's status 

182 shown in equation (8). Thus, what ecologists call throughflow Tj is a special case of Hubbell's status index 

183 Si when the model is deflned as in section (2.1.1). 
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184 2.2.1. Eigenvector and Environ Centrality 

185 To highlight its distinctiveness, Tj is contrasted with two alternative global centrality measures: eigen- 

186 vector centrality and environ centrality. As mentioned in the introduction, eigenvector centrality (EVC) 

187 describes the stable distribution of pathways, or when weighted as in flow networks the stable distribution 

188 of flow, passing through the nodes (Bonacich, 1972; Borgatti, 2005). In the context of directed flow networks, 

189 Fann and Borrett (2012) suggested using the average of the left w and right v hand eigenvector associated 

190 with the dominant eigenvalue of G to capture both the input and output, such that 



191 



EVC = [EVa] = i!^!i±^. (9) 



Note, in this calculation w and v are assumed to have been normalized so that their sum equals 1, which 

192 also implies that X]r=i ^^^i — 1- symmetric networks like those for which the eigenvector centrality was 

193 flrst defined Vi — Wi and averaging is not necessary. In directed flow networks Vi ^ Wi, and EVC captures 

194 the input and output oriented flows intensities. 

195 Fann and Borrett (2012) introduced average environ centrality (AEC) and argued that it is a better 

196 centrality indicator for ecosystem flow networks in part because it captures both the equilibrium dynamics 

197 (like EVC) and transient dynamics that occur along the initial shorter pathways in equation (5). This 

198 is important because in highly dissipative systems like trophic ecosystems, a large fraction of the total 

199 transactions might occur in these shorter pathways. Specifically, Borrett et al. (2010) found that in nine 

200 trophic ecosystem models 95% of TST required at most paths of length nine. AEC is defined as 



EC"' = [ecT] = -'"^ ' 



^=l 2Zj=i n 

En 
, _ i=l ""^ij 



Z^i=l Z^j = l ""1-3 

ECf + £:cr* 

AEC = [aecfc] (10) 

201 Although AEC is an improvement on EVC, both measures still suffer from two problems. The first 
issue is that the calculations required for EVC and AEC are not intuitive, which could be a barrier to their 
wider use in ecology (Fawcctt and Higginson, 2012). The second more substantive issue is that they fail to 
recognize or capture the external environmental forcing occurring in these open systems. Both measures 
are built on the non-dimensional fiow intensity matrices that represent the potential fiows or the flows if 

206 each node had a unit input. However, to recover the realized or observed system activity these matrices 

207 must be multiplied by the boundary vector as in equation 6 (see HubbcU, 1965). A critical issue is that 

208 the vector of boundary inputs in ecosystem models tends to be highly heterogeneous (Borrett and Freeze, 

209 2011), which differentially excites the potential flow pathways captured in G and N. Given these issues, in 

210 many applications Tj is a better indicator of the functional importance of a node because its calculation is 

211 more intuitive and because it captures the system's environmental forcing. 

212 The difference between these indicators can be substantive as illustrated for the Ythan Estuary and 

213 Chesapeake Bay ecosystem models (Fig. 3). In the Ythan Estuary, T is highly rank correlated with EVC 

214 (Spearman's p ~ 0.79) and AEC (p — 0.82), but T ranks the Nutrient Pool, Suspended POC, and Benthic 

215 Macrophytes as the top three nodes, which is not the case for the other two indicators. The flrst two of these 

216 nodes have boundary input. The Spearman rank correlation between T and EVC and AEC is generally 

217 less in the Chesapeake Bay model (p — 0.22 and p — 0.55, respectively). Again, EVC and AEC discount 

218 the importance of some nodes. In this case, the top three nodes - Phytoplankton, Suspended Particulate 

219 Carbon, and Dissolved Organic Carbon - have non-zero boundary input. Thus, T better captures the 

220 importance of nodes that connect the system to its external environment, and how this influence propagates 

221 throughout the system. 

222 In summary, throughflow is a global centrality indicator of the functional importance of nodes in a flow 

223 network. It is a special case of what Hubbell (1965) deflned as a status score in sociology. Due to the natural 
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224 discounting of longer pathways as energy or matter dissipates from the system, it has the desirable properties 

225 of mesoscale centrality measures advocated for by Jordan ct al. (2006). While it is similar to eigenvector 

226 and environ centrality measures, it is more intuitive to calculate and better captures the environmental 

227 forcing of the internal system activity. The next section applies T centrality to characterize the distribution 

228 of functional importance in 45 ecosystem models. 

229 3. Application — Materials and Methods 

230 Given that Tj is a global indicator of an ecosystem component's functional importance, we can now 

231 investigate the distribution of this importance in ecosystems. 

232 3.1. Ecosystem Model Database 

233 I applied flow analysis to 45 trophic ecosystem models selected from the literature and calculated Tj to 

234 investigate the throughflow centrality distributions (Table 1). To be included in this data set, the models 

235 needed to have at least 10 compartments, have a food web at their core (i.e., trophic models), and be 

236 empirically-based in the sense that the original modelers were attempting to represent a real ecosystem and 

237 used empirical measurements to parametrize part of the fluxes. If two models existed in the literature for 

238 the same system, only the least aggregated model (higher n) was included. Ten (22%) of these models are in- 

239 eluded in Dr. Ulanowicz network collection on his website (http://www.cbl.uniccs.edu/~ulan/ntwk/network.html). 

240 This data set also overlaps 80% with the models recently analyzed for resource homogenization (Borrctt and Salas, 

241 2010), dominance of indirect effects (Salas and Borrctt, 2011), and environ centrality (Fann and Borrett, 

242 2012). The full set of models are available at http://people.uncw.edu/borrctts/research.html. Forty-four 

243 percent of the models were not initially at steady-state, and were therefore balanced using the AVG2 algo- 

244 rithm (Allesina and Bondavalli, 2003). 

245 3.2. Centrality Comparison 

246 Rank correlation between T and AEC and EVC are shown for the Oyster Reef and Chesapeake Bay 

247 ecosystem models in section 2.2.1. Here, this result is generalized by examining distributions of the Spearman 

248 rank correlation between these measures in all 45 models in our database. 

249 3.3. Thresholds, and Dominants 

250 To characterize the T distributions within a model, I defined three thresholds. N50 is the number of 

251 nodes required to cumulatively account for 50% of TST when the compartments are rank ordered based on 

252 throughflow (largest to smallest). If a Monod function fit the cumulative fiow distribution, N50 would be 

253 equivalent to the half saturation constant. Ngo and N95 are the number of nodes required to recover 80% 

254 and 95% of TST. 

255 These thresholds are illustrated for the Bothnian Sea, Chesapeake Bay, and Sylt-R0m0 Bight ecosystems 

256 (Fig. 4). In the Bothnian Sea, only three nodes are required to generate 50% of the TST (N50 = 3), while 6 

257 and 8 nodes are required to account for 80% and 95% of TST, respectively (Ngo = 6 and Ngs =8). In the 

258 Chesapeake Bay model, these thresholds were N50 = 3, Ngo = 6, and N95 = 12, and in the Sylt-R0m0 Bight 

259 they were N50 — 3, Ngg = 7, and Ngs = 13. 

260 As the three models shown here have different numbers of compartments, n, it is difficult to compare these 

261 thresholds directly. For better comparisons, I normalized the thresholds by the model size as ^x/n * 100%. 

262 This gives the percent of nodes required to achieve the x% of TST. Fig. 4 shows that 33% of the model 

263 nodes are required to account for 95% of TST in the Chesapeake Bay model while only 22% of the nodes 

264 were required in the Sylt-R0m0 Bight model. This might be interpreted as indicating that system power is 

265 more concentrated in the Sylt-R0m0 Bight model. 

266 There are many ways of defining dominant species or compartments in ecological systems (e.g., Whittaker, 

267 1965; Fann and Borrett, 2012). Here, dominant compartments in the ecosystem were defined as the smallest 

268 subset of nodes required to recover 50% of TST. This definition lets us investigate both how many nodes are 

269 required for this (iVso) as well as their identity. For analysis, these compartments were classified as primary 
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270 producers (e.g., phytoplankton, submerged vegetation), dead organic matter (e.g., particulate organic mat- 

271 ter, dissolved organic matter), bacteria (e.g., free living bacteria, bacteria, bentliic bacteria), or other (e.g., 

272 filter feeders, meiofauna, large copepods). Detritus is technically a mixture of decomposers (some bacteria) 

273 with dead organic matter. For this analysis, detritus was grouped with the Dead Organic Matter. 

274 4. Results 

275 C'entrality Comparison 

276 As expected, EVC and AEC tend to be well correlated with T (Fig. 5). The median Spearman rank 

277 correlation between T and EVC is 0.69, with the values ranging between 0.11 and 0.87. Throughflow 

278 centrality is similarly correlated with AEC with a median value of 0.69. The distribution is visibly shifted 

279 to the right and has values ranging from 0.28 to 0.92. Notice that in no case is there 100% agreement or 

280 disagreement. 

281 4.2. Thresholds 

282 Figure 6 shows the cumulative flow development thresholds {N50, Ngo/n, Ng^/n) for the 45 trophic 

283 network models. There are several trends to note. First, the maximum number of nodes necessary to 

284 account for 50% of TST was 4. While in the Bothnian Bay ecosystem model this is 33% of the nodes, it 

285 is only 3.2% of the nodes in the Florida Bay model. Second, as the models increase in size (n) both Nso/n 

286 and Ng^/n tend to decline. Third, Figure 6b shows that in the majority (73%) of the models, 20% of the 

287 nodes or fewer account for 80% or more of the system activity. 

288 4-3- Dominants 

289 Figure 6a shows that 4 or fewer nodes are required to account for 50% of the TST and thus meet the 

290 criteria as dominants. The majority (46%) of the models analyzed had three dominant nodes, while another 

291 29% had only two dominant compartments (Fig. 7a). 

292 Table 2 identifies the 130 dominant nodes in each of the 45 ecosystems. The authors of the original 

293 models did not necessarily use identical categorizations for different ecosystem components, but it is possible 

294 to classify the compartments into four functional groups: primary producers, dead organic matter, bacteria, 

295 and a final category for anything else (other). Figure 7b shows the fraction of models that had at least one 

296 dominant in each of these categories. Thus, 82% of the models had at least one dominant compartment 

297 that functioned as a primary producer; 91% had a dominant compartment that was categorized as dead 

298 organic matter. Bacteria were also common. Only 9 of the dominant nodes did not fall into one of these 

299 three categories, and they only appeared in 7 of the models. 

300 5. Discussion 

301 Next I consider the theoretical development and its initial ecological application presented in this paper 

302 from three perspectives. First, I highlight some of the advantages and disadvantages of recognizing that 

303 system throughflow is a centrality indicator. Second, I contemplate the import of this discovery for un- 

304 derstanding ecological system organization, growth, and development. Third, I identify additional possible 

305 applications of this innovation. 

306 5.1. Throughflow as a Centrality 

307 A primary contribution of this paper is to recognize that throughflow T, a measure used by ecologists 

308 for some time (e.g., Finn, 1976; Ulanowicz, 1986; Fath and Patten, 1999), is a centrality measure as defined 

309 in the social science (Hubbell, 1965; Fricdkin, 1991; Wasscrman and Faust, 1994) and now used in general 

310 network science (Brandos and Erlebach, 2005). An advantage of connecting throughflow and centrality is 

311 that ecologists can now access, apply, and further develop the existing body of work on centrality. For 

312 example, many centrality measures have been proposed, but sociologists can generally classify them into 
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313 one of three types (Freeman, 1979; Friedkin, 1991; Wasserman and Faust, 1994; Borgatti and Everett, 2006). 

314 The first type are degree based measures. These measures can vary in the size of the neighborhood considered 

315 - from the immediate local neighborhood to global measures that consider the whole system (e.g., Estrada, 

316 2010). This type of centrality is generally interpreted as the influence of the node on the network activity 

317 or its power to change the activity (Bonacich, 1987). A second type of centrality is termed closeness and 

318 is based on the shortest paths or geodesic distances between nodes. Friedkin (1991) suggests that these 

319 measures indicate the immediacy of a node's ability to influence the network. A third commonly described 

320 type of centrality is betweeness (Freeman, 1979; Freeman et al., 1991). A node's betweeness centrality is 

321 its importance in transmitting activity between individuals or subgroups in the network. Thus, there is a 

322 recognition of several different but complementary ways in which individuals in a system can be central. 

323 In this broader context of centralities, Hubbell's status is a global, weighted, degree based centrality that 

324 is typically interpreted as the node's influence on the whole system activity or its power to change the whole 

325 system activity (Borgatti, 2005; Brandes and Erlebach, 2005). The formulation allows the node's centrality 

326 to be recursively changed by the centrality of the other nodes in the system as its walk connectivity is 

327 extended. Although Hubbell (1965) initially considered a potentially heterogeneous set of exogenous inputs, 

328 in practice a uniform set of inputs are typically used to consider the potential centrality. This is similar 

329 to the "unit" input analytical approach often used in network environ analysis (Fath and Patten, 1999; 

330 Whipple et al., 2007; Borrctt and Freeze, 2011). In the ecological application of Hubbell centrality, the 

331 realized throughflow centrality is obtained using the observed exogenous inputs. 

332 Ecologists can further benefit from the sociologists previous applications of centrality. For example, 

333 Hubbell initially used his centrality as a tool to detect subcommunities or cliques within the system. As 

334 this is again a common concern for ecologists (Pimm and Lawton, 1980; Allesina et al., 2005; Borrctt ct al., 

335 2007), we may be able to utilize his procedure to address this problem in the future. This would follow 

336 Krausc et al.'s (2003) successful application of a different social network analysis clique finding algorithm 

337 to food webs. 

338 Another advantage is that we may be able to recognize other ENA measures as centrality type indicators. 

339 For example, several of Friedkin's (1991) descriptions of alternative centrality measures for what he called 

340 "total effects centrality" were very similar to what Whipple et al. (2007) called total environ throughflow 

341 (TET). Thus, TET may also be a type of weighted degree centrality measure that indicates the relative 

342 contribution of each environ to the whole system activity. Hincs et al. (2012) has already begun to explore 

343 this possibility while investigating nitrogen cycling model of the Cape Fear River estuary 

344 There are two potential disadvantages of recognizing throughflow as a centrality indicator. First, it 

345 could contribute to the proliferation of centrality measures that can be overwhelming. This has led to 

346 multiple papers trying to identify the unique contributions of specific indicators amongst a set of competing 

347 indicators (e.g., Newman, 2006; Jordan et al., 2007; Valente et al., 2008; Bauer ct al., 2010; Baranyi ct al., 

348 2011). In this case, however, I argue that we are not creating a new centrality index to add to the confusion, 

349 but identifying that a commonly calculated measure is a form of an existing centrality measure. A second 

350 disadvantage might be that the current use and implementation of Hubbell's centrality available in software 

351 packages may be simplified from its original formulation, as appears to be the case in Ucinet (Borgatti et al., 

352 2002). The output of the Hubbell centrality analysis in Ucinet does not match the throughflow vector as 

353 calculated with NEA.m (Fath and Borrett, 2006) 

354 As expected, T generally correlates well with average eigenvector centrality (EVC) and average environ 

355 centrality (AEC) for the 45 models examined. This suggests that these different global degree-based central- 

356 ity measures capture some of the same information about the relative importance of the nodes for the system 

357 function. However, the correlations were variable - in some cases the rankings were quite different (e.g., 

358 median Spearman correlations were 0.69 and the lowest was 0.11) suggesting that each measure captured 

359 some unique information. Examining both the formulation of the three centrality measures as well as the 

360 example in Figure 3, a key difference is that T captures the importance of a node for connecting the system 

361 to the external world. For example in the Ythan estuary model, the Nutrient Pool and Suspended POC 

362 both have large inputs that contribute to their importance in T. Thus in applications where the boundary 

363 inputs are an important consideration, an indicator like throughflow centrality may be the best choice. For 

364 example, Borrett and Freeze (2011) argued that this system-environment coupling is critical for ecologists 
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365 and environmental scientists even when the analytical focus is on the within system environments. 

366 5.2. Throughflow and Ecosystem Organization and Development 

367 Ecologists have a long interest in the organization, growth, and development of ecosystems (e.g., Odum, 

368 1969; Ulanowicz, 1986; J0rgensen et al., 2000; Gundcrson and Holling, 2002; Lorcau, 2010). What are 

369 the processes that create, constrain, and sustain ecological systems? Scientists investigating this prob- 

370 lem have hypothesized a number of goal functions or orientors that might guide the growth and develop- 

371 ment of these self-organizing systems (Schneider and Kay, 1994; Miiller and Leupelt, 1998; J0rgensen et al., 

372 2007). Hypothesized orientors include the tendency for ecosystems to maximize power (Lotka, 1922; 

373 Odum and Pinkerton, 1955), maximize biomass or storage (J0rgensen and Mejer, 1979), maximize dissipa- 

374 tion (Schneider and Kay, 1994), and maximize emergy (Odum, 1988). Fath et al. (2001) used the network 

375 framework to show how these different orientors can be complementary. 

376 Patten (1995) suggested that throughflow in network models of energy flux can be interpreted as a 

377 measure of power in a thermodynamic sense. He argued that TST indicates the total power output of an 

378 ecological system. This operationalized Lotka's (1922) maximum power principle for evolutionary systems 

379 and Odum and Pinkcrton's (1955) hypothesis that ecological systems tend to maximize their power in a 

380 network context. Given this interpretation of TST, Tj is therefore the partial power of each node (j = 

381 1, 2, . . . , ri) in the network. Interestingly, this thermodynamic interpretation to throughflow aligns with the 

382 social interpretation of this type of centrality as the power to influence the system (Bonacich, 1987). 

383 Recognizing that network nodes in ecosystem models represent subsystems in a hierarchical context 

384 (Allen and Star, 1982), then we can extend the maximum power hypothesis to each node. As all nodes 

385 would experience the same attraction to increase Tj, we might expect the TjS to be more similar (towards a 

386 uniform distribution). However, this maximization remains restrained by the evolutionary constraints of the 

387 individual organisms, including their participation within the existing ecosystem (Walsh and Blows, 2009; 

388 Guimaracs Jr et al., 2011). For example, Ulanowicz (1997, 2009) argues that the formation of autocatalytic 

389 cycles can be an agency for ecosystem growth and development. These cycles can provide the positive 

390 feedback and selective pressure for individual nodes to tend to increase their Tj. They also provide a 

391 selection pressure such that alternative nodes within an autocatalytic cycle compete for participation in 

392 throughflow and can be replaced by higher performing entities. Ulanowicz (1997) further argues that the 

393 tendency of these cycles for centripitality - in this context attracting and capturing more resources - leads to 

394 the emergence of a system autonomy from the material cause of the system. Thus, evolutionary constraints 

395 on species and the system constraints of interacting autocatalytic cycles might increase the variability of Tj 

396 despite the homogenizing effect of the tendency to maximize throughflow. 

397 The throughflow threshold analysis of the 45 ecosystem models presented here indicates that throughflow 

398 centrality is far from uniform as it appears to follow something more like Pareto's 80-20 rule in which 80% 

399 of the activity is done by 20% of the group (Reed, 2001). This suggests that throughflow centrality may 

400 be similar to if not exactly the scale free degree distributions commonly found in other types of complex 

401 systems (Barabasi, 2002). In addition, all but 8 of the dominant or most central nodes could be classified as 

402 primary producers, dead organic matter, or bacteria. This aligns with what we might expect from ecosystem 

403 theory in general and the importance of autocatalytic hypercycles like the autotroph -H- decomposer cycle 

404 (Ulanowicz, 1997; Wilkinson, 2006). 

405 5.3. Applications 

406 Network modeling and analysis, Input-Output Analysis, and material flow analysis have broad applica- 

407 tion. The ideas originated in macro economics (Lcontief, 1966) and as has been discussed are used in both 

408 sociology and ecology. Thus, throughflow centrality may be useful in multiple domains of inquiry. 

409 Beyond the theoretical considerations for ecosystem growth and development, there are a number of 

410 ways in which the throughflow centrality indicator could be usefully applied for ecosystem management, 

411 conservation, and restoration. For example, the throughflow centrality analysis suggests which species or 

412 groups of species should be targeted in the goal is to increase or decrease the system activity. The impact 

413 of manipulating a more central node should be greater than modifying a less central node. 
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414 Materials flow analysis is an important tool for industrial ecology (Bailey et al., 2004a,b; Suh and Kagawa, 

415 2005; Gondkar ct al., 2012) and urban metabolism (Kennedy et al., 2011; Zhang ct al., 2012; Chen and Chen, 

416 2012). The specific ENA methods described in this paper have been used to analyze the sustainability of ur- 

417 ban metabolisms (Bodini and Bondavalli, 2002; Zhang et al., 2010; Chen and Chen, 2012). Chen and Chen 

418 (2012) shows how throughflow can be grouped according to compartment "trophic levels" to build produc- 

419 tivity pyramids for cities that are then comparable to expected trophic productivity pyramids in ecology. 

420 Thus, the recognition that T is a centrality indicator could have a broad utility for these disciplines. 

421 ENA is an ecoinformatic tool and shares many goals and characteristics with network analysis in the 

422 field of Systems Biology. For example, Hahn and Kern (2005) showed that genes with higher centrality 

423 tend to be functionally more important in protein-protein interaction networks. While thermodynamically 

424 conserved flows are not normally the focus of the systems biology network models (omics) making it difficult 

425 to apply the flow analysis and ENA more broadly, Kritz et al. (2010) suggest a way of liking a metabolic 

426 network model to the underlying chemical fluxes and reactions. If this technique proves robust, then the 

427 throughflow centrality might be useful in this domain as well. 

428 6. Conclusions 

429 In summary, this paper makes two primary contributions. First, I show that throughflow (Tj) in network 

430 input-output models is a global indicator of the relative importance or power of each node in the network 

431 with respect the whole system activity. As calculated in ecological network analysis, this is a special case 

432 of Hubbell centrality (Hubbell, 1965). Second, when applied to trophic network models of ecosystems, 

433 throughflow centrality shows the tendency of this power to be concentrated in a small set of nodes that tend 

434 to categorized as primary producers, dead organic material, or bacteria. This is consistent with previous 

435 theory regarding the growth and development of ecological systems. 

436 To address the wicked problems (Rittel and Webber, 1973) of our time like economic challenges and 

437 global climate change, we will need to be both creative and innovative. An innovation in this paper is to 

438 join the throughflow concept in flow analysis and the centrality concept developed in the social sciences. I 

439 expect this to be a useful union that will enable new analysis and management of complex systems of many 

440 kinds including urban metabolisms, industrial ecosystems, and biogeochemical cycling and trophic dynamics 

441 in natural ecosystems. 
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Table 1: Forty-five empirically-based trophie ecosystem network models. 



Model 


units 


nt 


ct 


Boundary'^ 


TSTt 




Source 


Bothnian Bay 


gC m yr 


12 


0.22 


44 


184 


0.23 


Sandbcrg et al. (2000) 


Bothnia!! Sea 


gC m yr 


12 


0.24 


117 


562 


0.31 


Sandberg et al. (2000) 


Ythan Estuary 


gC m yr 


13 


0.23 


1,259 


4,182 


0.24 


Baird and Milne (1981) 


O J 1 TV T f ' ' \ 

bundarban Mangrove (virgm) 


1 1 —2 —1 
kcal m ^ yr ^ 


14 


0.22 


117,959 


441,214 


0.16 


Ray (2008) 


O 11 TVT /I" 1\ 

bundarban Mangrove (reclaimed) 


1 1 —2 —1 
kcal m yr 


14 


0.22 


38,485 


103,057 


0.05 


Ray (2008) 


Baltic Sea 


—2 J — 1 

mgC m d 


15 


0.17 


603 


1,974 


0.13 


Baird et al. (1991) 


Ems Estuary 


/~i —2 J — 1 
mgC m d 


15 


0.19 


283 


1,067 


0.32 


Baird et al. (1991) 


O J 1 1 T T 11" 

Southern Benguela Upwelling 


1-* —2 J — 1 
mgC m d 


16 


0.23 


715 


2,546 


0.31 


Baird et al. (1991) 


Peruvian Upwelling 


mgC m ^ d ^ 


16 


0.22 


14,928 


33,491 


0.04 


Baird et al. (1991) 


Crystal River (control) 


—9 1 — 1 

mgC m d 


21 


0.19 


7,358 


15,063 


0.07 


Ulanowicz (1986) 


Crystal River (thermal) 


—2 J — 1 

mgC m d 


21 


0.14 


6,018 


12,032 


0.09 


Ulanowicz (1986) 


JUT 1 T 

Charca de Maspalomas Lagoon 


—2 J — 1 
mgC m d 


21 


0.12 


1,486,230 


6,010,331 


0.18 


Almuma et al. (1999) 


Northern Benguela Upwelling 


r~< —2 J — 1 

mgC m d 


24 


0.21 


2,282 


6,609 


0.05 


Heymans and Baird (2000) 


Swartkops Estuary 


—2 J — 1 

mgC m ^ d ^ 


25 


0.17 


2,860 


8,950 


0.27 


nil 1 ■ 1 /<~\nr\r'\ 

bcharler and Baird (2005) 


Sundays Estuary 


—0 1 — 1 

mgC m d 


25 


0.16 


4,442 


11,940 


0.22 


Scharler and Baird (2005) 


Kromme Estuary 


mgC m d 


25 


0.16 


2,571 


11,088 


0.38 


bcharler and Baird (2005) 


Neuse tstuary (early summer 1997) 


—2 J — 1 
mgC m d 


30 


0.09 


4,385 


13,828 


0.12 


Baird et al. (i004b) 


Neuse tstuary (late summer 1997) 


—2 J — 1 

mgC m d 


30 


0.11 


4,640 


13,036 


0.13 


Baird et al. (2004b) 


Neuse tstuary (early summer 1998) 


—2 A — 1 
mgC m ^ d 


30 


0.09 


4,569 


14,025 


0.12 


Baird et al. (i004b) 


Neuse tstuary (late summer 1998) 


—2 J — 1 
mgC m d 


30 


0.1 


5,641 


15,032 


0.11 


Baird et al. (i004b) 


Gulf of Maine 


— 2 —1 
g WW m yr 


31 


0.35 


5,054 


18,382 


0.15 


Link et al. (2008) 


Georges Bank 


— 2 —1 

g WW m yr 


31 


0.35 


4,381 


16,890 


0.18 


Link et al. (2008) 


Middle Atlantic Bight 


— 2 —1 

g WW m yr 


32 


0.37 


4,869 


17,917 


0.18 


Link et al. (2008) 


Narragansett Bay 


r-i 9 1 

mgC m ^ yr ^ 


32 


0.15 


693,846 


3,917,246 


0.51 


Monaco and Ulanowicz (1997) 


Southern New England Bight 


— 2 — ! 
g WW m yr 


QQ 
OO 


U.oo 




1 7 C^OV 

1 /,oy / 


n 1 ^1 
U. ib 


LinK et al. (zuuo ) 


Chesapeake Bay 


mgC m~2 yr~^ 


36 


0.09 


888,791 


3,227,453 


0.19 


Baird and Ulanowicz (1989) 


St. Marks Seagrass, site 1 (Jan.) 


mgC m~2 d~^ 


51 


0.08 


515 


1,316 


0.13 


Baird et al. (1998) 


St. Marks Seagrass, site 1 (Feb.) 


mgC m~2 d~^ 


51 


0.08 


602 


1,591 


0.11 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Jan.) 


mgC m~2 d~^ 


51 


0.07 


603 


1,383 


0.09 


Baird et al. (1998) 


St. Marks Seagrass, site 2 (Feb.) 


mgC m-2 d-l 


51 


0.08 


801 


1,921 


0.08 


Baird et al. (1998) 


St. Marks Seagrass, site 3 (Jan.) 


mgC m-2 d-l 


51 


0.05 


7,809 


12,651 


0.01 


Baird et al. (1998) 


St. Marks Seagrass, site 4 (Feb.) 


mgC m-2 d-^ 


51 


0.08 


1,433 


2,865 


0.04 


Baird et al. (1998) 


Sylt R0m0 Bight 


mgC m-2 d-^ 


59 


0.08 


683,448 


1,781,029 


0.09 


Baird et al. (2004a) 


Graminoids (wet) 


gC m— 2 yr— ^ 


66 


0.18 


6,272 


13,677 


0.02 


Ulanowicz et al. (2000) 


Graminoids (dry) 


gC m— 2 yr— ^ 


66 


0.18 


3,473 


7,520 


0.04 


Ulanowicz et al. (2000) 


Cypress (wet) 


gC m— 2 yr— ^ 


68 


0.12 


1,419 


2,572 


0.04 


Ulanowicz et al. (1997) 


Cypress (dry) 


gC m— 2 yr— ^ 


68 


0.12 


1,036 


1,919 


0.04 


Ulanowicz et al. (1997) 


Lake Oneida (pre-ZM) 


gC m— 2 yr— ^ 


74 


0.22 


1,035 


1,698 


0.00 


Miehls et al. (2009a) 


Lake Quinte (pre-ZM) 


gC m— 2 yr— ^ 


74 


0.21 


989 


1,518 


0.00 


Miehls et al. (2009b) 


Lake Oneida (post-ZM) 


gC m— 2 yr" ^ 


76 


0.22 


811 


1,463 


0.00 


Miehls et al. (2009a) 


Lake Quinte (post-ZM) 


gC m— 2 yr— ^ 


80 


0.21 


1,163 


2,108 


0.01 


Miehls et al. (2009b) 


Mangroves (wet) 


gC m— 2 yr— ^ 


94 


0.15 


1,532 


3,266 


0.10 


Ulanowicz et al. (1999) 


Mangroves (dry) 


gC m— 2 yr— ^ 


94 


0.15 


1,531 


3,272 


0.10 


Ulanowicz et al. (1999) 


Florida Bay (wet) 


mgC m— 2 yr— ^ 


125 


0.12 


739 


2,721 


0.14 


Ulanowicz et al. (1998) 


Florida Bay (dry) 


mgC m— 2 yr" ^ 


125 


0.13 


548 


1,779 


0.08 


Ulanowicz et al. (1998) 



T n is the number of nodes in the network mode the model connectance when L is the number of direct links or 

energy— matter transfers, TST — "^^2 fij + -^i is the total system throughftow, and FCI is the Finn Cycling Index (Finn, 1980). 
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Table 2: Dominant ecosystem components as identified by throughflow centrality with primary producers labeled with a green 
box, colored in a brown box with white letters, and Jjacteria in a pink box. These are the model nodes 



required to generate 50% of total system throughflow (N50). 



T2 



T3 
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an Sea 






Vlacrofauna 




Ythan 


Estuary 








Suspended POC 1 


Sunda 
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virgin 






Macrophytes 
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rcelai 
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Detritus 
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Baltic 


Sea 






Pelagic Production"^ 


Mcsozooplankton 



Ems Estuary 

Soutliorn Bengucla Upwelling 
Peruvian Upwelling 
Crystal River (control) 
Crystal River (thermal) 
Charca de Maspalomas Lagoon 
Northern Bengucla Upwelling 
Swartkops Estuary 
Sundays Estuary 
Kromme Estuary 

Ncuse Estuary (early summer 1997) 
Neuse Estuary (late summer 1997) 
Ncuse Estuary (early summer 1998) 
Ncuse Estuary (late summer 1998) 



■■l[U.lj.l.lJ.»J.IJI 
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Macrophytes 

Macrophytes 
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mm. 



1 
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Everglade Graminoids (dry) 

Cypress (wet) 

Cypress (dry) 
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Lake Oneida (post-ZM) 

Lake Quinte (post-ZM) 

Mangroves (wet) 

Mangroves (dry) 

Florida Bay (wot) 
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Diatoms 
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Figure 1: Hypothetical network model (a) with its associated (b) degree and (c) eigenvector centrality. Degree centrality is a 
local measure while eigenvector centrality is a global indicator of node importance. 
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Figure 2: Rank order throughflow centrality for the Gulf of Main ecosystem. 
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Figure 3: Comparison of throughfiowr centrality (TC), average eigenvector centrality (EVC), and average environ centrality 
(AEC) in the Ythan Estuary (a) and Chesapeake Bay (b) ecosystem networks. Model compartments are rank ordered by TC. 
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Figure 4: Rank ordered throughflow (a, b, c) and cumulative throughfiow (d, e, f) for the Bothnian Sea (a,d), Chesapeake 
Bay (b,e), and Sylt-R0m0 Bight network models (c,f). Throughflow has the units shown in Table 1. The thick horizontal line 
in a, b, and c shows what throughflow would be if each node contributed equally. The vertical lines in (c), (d), and (f) show 
the nodes at which 50% (N50), 80% (Ngo), and 95% (N95) of the total system throughflow is achieved. For the Bothnian Sea, 
these thresholds are achieved at 3, 6, and 8, respectively. In the Chesapeake Bay they are 3, 6, and 12, while in the Sylt-R0m0 
Bight they are 3, 7, and 13. 
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(a) Eigenvector Centrality 
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Figure 5: Histogram of Spearman rank correlation coefficients between throughflow centrality Tj and (a) average eigenvector 
centrality (EVC) and (b) average environ centrality (AEC) in 45 ecosystem models. 
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Figure 6: Rank order cumulative throughflow thresholds in 44 empirically based ecosystem models (models ordered by n with 
smallest at the top): (a) number of nodes required to account for 50% (N50), (b) percent of model nodes required to achieve 
80% (Ngo/n * 100%), and (c) 95% (Ngs/ra * 100%) of total system throughflow. 
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Figure 7: Dominant analysis: (a) the frequency of the 45 models with 1, 2, 3, or 4 dominant nodes (iVso), and (b) the percent 
of models with at least one dominant node in the three functional categories. 
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